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Abstract 

We propose an experiment to search for invisible decays of orthopositronium (o-Ps) with a 90% 
confidence sensitivity in the branching ratio as low as 10 -8 . Evidence for this decay mode would 
unambigously signal new physics: either the existence of extra-dimensions or fractionally charged 
particles or new light gauge bosons. The experimental approach and the detector components of the 
proposed experiment are described. 
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1 Based on a talk given at Workshop on Positronium Physics, Zurich, Switzerland, 30-31 May 2003. 



1 Motivation 



The Standard Model (SM) is very succesful, but is not yet a complete theory. There are several fun- 
damental questions, which require physics beyond the SM for a solution. Among them are e.g. the 
hierarchy problem, the origin of the charge quantisation and the dark matter composition. Although the 
new physics is expected typically at high energies, some models predict phenomena at low energies that 
might be observed in rare decays of the positronium, e.g. into a photon and an axion, a new particle 
introduced for the solution of the strong CP problem, for a review see[D[2]- 

Positronium, the positron-electron bound state, is the lightest known atom, which is bound and 
self-annihilates through the same, electromagnetic interaction. At the current level of experimental and 
theoretical precision this is the only interaction present in this system, see e.g. This feature has 

made positronium an ideal system for testing the accuracy of the QED calculations for bound states, in 
particular for the triplet (l 3 Si) state of Ps, orthopositronium (o — Ps). Due to the odd-parity under C- 
transformation o — Ps decays predominantly into three photons. As compared to the singlet (1 So) state 
(parapositronium), the long lifetime of o — Ps decay, (due to the phase-space and additional a suppression 
factors) gives an enhancement factor ~ 10 3 , in the sensitivity to an admixture of new interactions, which 
are not accommodated in the SM. 

This paper is focused on the search for o — Ps — > invisible decay, i.e. a process which is not accom- 
panied by energy deposition in a detector. Within the SM orthopositronium can decay invisibly into a 
neutrino-antineutrino pair. The o — Ps — > v e v e decay occurs through W exchange in the i-channel and 
e + e~ annihilation via Z . The decay width is [3] 

T(o - Ps -> v e v e ) « 6.2 • l(T 18 r(o -Ps^ 37) (1) 

For other neutrino flavours only the Z-diagram contributes. For I ^ e the decay width is |3] 

r(o - Ps -> vm) w 9.5 • icr 21 r(o - Ps 37) (2) 

Thus, in the SM the o — Ps — ► vv decay rate is very small and evidence for invisible decays would unam- 
biguously signal the presence of new physics. It may be worthwhile to remember that the process with 
analogous experimental signature, Z — > invisible decay plays a fundamental role in the determination of 
the number of lepton families. 

The new models that are relevant to the o — Ps —> invisible decay mode predict the existence either of 
i) extra-dimensions or ii) fractionally charged particles 5 , or iii) a new light vector gauge bosonpQ, 

or iv) dark matter of the mirror matter type|HllZj. To test these models, the required sensitivity in the 
branching ratio Br(o — Ps — » invisible) has to be at least as low as 10~ 8 . Note however that the search 
for the mirror matter requires in addition the experiment to be performed in vacuum (for more details 
see 0). Such an experiment is discussed and proposed in pT)j . 

The first experiment to search for invisible decay channels of o-Ps was performed a long time ago [B]- 
The best present limit for the branching ratio, set in 1992 by Mitsui et al.JS] is 

Br(o -Ps^ invisible) < 2.8 x 1(T 6 (3) 

The paper is organized as follows. In the next sections the physics of items i) — iii) is briefly discussed. 
The experiment, its setup components and the expected sensitivity are discussed in Sections 2, 3 and 4 
respectively. Section 5 contains concluding remarks. 

1.1 Extra dimensions 

Models with infinite additional dimensions of the Randall-Sundrum type (brane-world models) with a big 

compactification radius ^5 [HI could provide the natural solution to the gauge hierarchy problem. 

Recently Dubovskj, Rubakov and Tinyakov 0] pointed out that in R-S brane world models with 
localized bosons, massive gauge bosons localized to the Plank brane are always unstable. The reason is 
that bulk modes can have arbitrarily small masses. Any massive mode localized to the brane is then 
kinematically allowed to decay into these modes. 
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In a recent paper [2] this mechanism was applied to calculate the decay of o-Ps into invisible bulk modes 
and obtained 

Br(o — Ps — > 7* — > additional dimension(s)) = (4) 

4(tt 2 -9) ' a 2 ' 16 1 fc ' ~ 1 fc J 

For a solution of the gauge hierarchy problem the parameter fc is expected to be k < O(10) TeV. Indirect 
measurements at LEP of the T(Z — > invisible) decay width constrain the parameter to fc > 2.7 TeV. 
Using this bound one can find 

Brio — Ps — > additional dimension(sj) < 0.4 • 10~ 9 (5) 

The direct LEP measurements of T(Z — > invisible) results in a less stringent limit 012] ■ 

i?r(o — Ps — > additional dimension(s)) < 8 • 10 -8 (6) 

These estimates show that the region of ~ 10~ 8 — 10 -9 for the branching ratio is of great interest for the 
observation of extra-dimensions. 



1.2 Millicharged particles 

In 1986 Holdom jS] showed that grand unified models can be constructed in a natural way, adding 
a second, unobserved, photon to the interaction Lagrangian. In this type of model particles with an 
electric charge very small compared to the electron, are predicted. If such milli-charged particles exist 
with a small mass, the o — Ps could decay apparently invisibly, since the particles would most likely 
penetrate any type of calorimeter without interaction. The corresponding decay width is JS] 



T(o-P S ^XX) = ^f^.k.F(^f), (7) 



where Qx is the electric charge of the X-particle (Q e = 1), k = 1, F(x) = (1 — ^x)(l — x) 2 for spin \ 
and k — j, F(x) = (1 — x)i for a spin-less X-particle. For spin | millicharged X-particles and nix *C m e 
one can find from the experimental bound (3) that Qx < 8.6 • 10 -5 [H]-A search for millicharged particles 
through the o — Ps — > invisible decay with the sensitivity of B>r(p — Ps — > invisible) ~ 10 -8 would touch 
the parameter space which is not excluded by results of the recent direct experiment at SLAC |16j . 

1.3 New light X-boson 

A light vector A-boson with the interaction Lagrangian 

L x = gxi>i^x» (8) 

will also lead to invisible decays of o — Ps. 1 Supposing that the A-boson interacts with other particles 
(fermions) or (as a consequence of the Higgs mechanism) with itself and a scalar the contribution of the 
A-boson to the electron anomalous magnetic moment is given by the well known formula 

5ae = 2*[ l 2 (9) 

T JO X 2 + (l- X )^ 

For mx <S m e from the bound arising from the determination of the fine structure constant with the 
Quantum Hall Effect (QHE)pQ, one can find that ax < 3 x 10~ 10 . For the case of the heavy A-boson 
(mx S> m e ) the bound on the anomalous electron magnetic moment leads to the bound 

2 

ax^r- <4.5-10~ 10 (10) 
m x 



1 For the recent phenomenological bounds in models with a light vector X-boson related to the muon [g — 2) and so-called 
NuTeV anomalies see, respectively El, [TS| and [15]. 
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For the reaction o — Ps — ► X* — > X\Xi (here X* is a virtual X-boson and -X^ is a fermion (sterile 
neutrino) or a scalar particle) one finds that 

Br(o - Ps -> X* -> = (11) 



4(7r 2 - 9) m; 



fc-F( — ^)(1- ., 

mri Z, fit 



O-Ps 



where F(x) has been defined before and axxi — 9xx 1^- F rom the bound of the QHE we find for 
rrix <S m e that 

Br(o -Ps^ X* -> XiXi) < k x 2 • 10~ 3 • a X x t (12) 

To have an experimentally interesting branching ratio of the order O(10~ 8 ), setting ax = 3 x 10~ 10 we 
must have axxi ~ 5 • 10 -6 . 

In the opposite limit mx S> m e the corresponding bound reads 

— TTl^ 

Brio - Ps -> X* -> X^i) < fc x 3 • 10~ 3 • a XXl ■ —§- (13) 

m x 

2 

For an X-boson mass close to the orthopositronium mass, we have the enhancement factor (1 — m "' x )~ 2 
in formula (11) for the branching ratio, and hence, the coupling constant ax • olxxi could be smaller. 



2 Proposed setup and experimental technique 

In this section the design of the experiment to search for o — Ps — > invisible with a sensitivity better than 
10~ 8 is presented. The main components of the detector (shown in Figs. 1, 2) are: the positron source 
( 22 Na), the positron tagging system, composed of a scintillating fiber viewed by two photomultiplicrs 
(PM), the positronium formation Si02 target and a hermetic 7-detector. The coincidence of the PM 
signals (see FigEJ from the positrons crossing the fiber, opens the gate for the data acquisition (DAQ). 
In the off-line analysis the 1.27 MeV photon, which is emitted from the source simultaneously with the 
positron, is required to be in the trigger BGO counter (shown in Fig[2J) resulting in a high confidence level 
of positron appearence in the positronium formation region. A positron, which enters the Si02 target 
may capture an electron creating positronium. The calorimeter detects, either the direct 27 annihilation 
in flight or the 2(3) photons from the para (ortho)-positronium decays in the target. 

The occurrence of the o — Ps — » invisible decay would appear as an excess of events with zero-energy 
deposition in the calorimeter above those expected from the Monte Carlo prediction (see section 4) or 
from the direct background measurement. 

This measurement presents a new feature of this type of experiment. The idea is to obtain a pure 
o-Ps decay energy spectrum by comparing two different spectra from the same target filled either with 
N2 (low o- Ps quenching) or with air, where the presence of paramagnetic O2 will quench the fraction of 
o-Ps in the target from 10% down to 3%, due to the spin exchange mechanism 20 : 

o-Ps + 0{]])^p-Ps + 0{U) (14) 

Thus the subtraction of these properly normalized spectra will result in a pure o-Ps annihilation energy 
spectrum in the 7-detector. 

In addition, compared to the previous experiments the region around the target has been de- 
signed with as little dead material as possible in order to reduce the photon absorption. The confidence 
level for tagging the positron appearance in the target and the efficiency of annihilation photon detection 
have been improved. 



2.1 Source of positrons 

The radioactive isotope 22 Na, which has a half-life of 2.6 years, was chosen as a source of positrons. The 
positron emission from the source is accompanied by the prompt emission of a 1.27 MeV photon (~ 3 
ps delay) (see Fig. 3). This signature provides an inefficiency for positron tagging less than ~ 10 -9 . The 
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Figure 2: Side view of the calorimeter. 



source activity was selected based on a compromise between required trigger rate and a decrease in the 
signal efficiency due to the overlap of events close in time (pile up effect). The optimal source activity has 
been determinated to be about 50 kBq, for which, according to the simulation, the signal efficiency drops 
down to 84% (the signal of o — Ps — > invisible is defined as an event depositing less than 50 keV in the 7- 
detector). In 10% of the cases the 1.27 MeV photon from the source is emitted due to the electron capture 
(EC), without positron emission. As discussed in the next section, this process introduces background 
related to the source container material, which was eliminated by depositing the radioactive isotope 
directly on the scintillating fiber. A comparison between 22 Na and other (3 + radioactive isotopes with 
a sufficiently long lifetime has been made to combine the required confidence level for positron tagging 
with the minimal amount of dead target material and a high stopping efficiency of positrons in the target. 
Taking all these considerations into account, we conclude that the 22 Na is the best compromise for our 
measurements. 



90% Positron with end point energy 
543 keV, 10% electron capture, 
no positron emission 



100% Photons 
energy 1275 keV 



- 2 Ne 




:: Na 



0.06% Positron with end 
point energy 1 890 ke V 



Figure 3: Decay scheme of 22 Na 
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2.2 The positron tagging system 



The tagging system for the experiment was designed in order to have a confidence level for the positron 
appearance better than 10~ 9 . It is composed of a thin scintillating fiber viewed by two PM's. It provides 
a relatively high efficiency and a small amount of dead material. 

2.2.1 The scintillating fiber 

For a good positron tagging efficiency , the attenuation length and the light yield of the fiber are crucial 
parameters. After testing the response of different scintillating fibers, the Bicron BFC 12 with a diameter 
of 0.5 mm has been selected. They have been prefered over the I mm fibers used in the experiment^ I j 
because they have less material and the possibility to have a more hermetic calorimeter, due to the smaller 
gaps between the crystals, through which the fiber has to penetrate to transport the signal outside the 
calorimeter volume. Similar to experiment |21| the fiber is squeezed in the middle to get a thickness 
of about 100 /im. The 22 Na source is deposited in the middle of the squeezed region. The measured 
efficiency for positron detection is slightly less than 50%. 

2.2.2 The photomultipliers 

The PM's XP 2020 have been chosen for their low noise and for their spectral sensitivity, that fits the 
emission spectrum of the fiber. The use of two PM signals in coincidence decreases the probability of 
accidental coincidences with the 1.27 MeV photon to a level of 3.2 x 10~ n (see section 3). In the future 
we plan to test Hamamatsu H3165-11 PM's with base R647, because of their high efficiency for single 
photon counting. 2 

2.3 Positronium production target 

A S1O2 aerogel hemisphere with a radius of 5 mm and a density of 0.1 (//cm 3 , where positronium can be 
formed with about 5% probability (as measured directly from the lifetime spectrum), is used as a target. 
The aerogel is a porous material, where the o-Ps migrating in the intergranular space can decay almost 
freely. The dimension has been optimised to minimize the dead material, while keeping, enough material 
to stop positrons efficiently. The simulation shows that the probability that more than 900 keV photons 
energy is absorbed in it, for the given dimension, is less than 10~ 10 . Before its installation the aerogel is 
evacuated in a vacuum chamber (10 -3 Torr) and backed out for 2 hours to remove the water from the 
pores. Then the target is filled cither with high purity Nitrogen or with air. 

2.4 Calorimeter 

The optimal choice of the 7-detector (ECAL) can be made by the following considerations. The total 
ECAL mass W is given roughly by 

W ~ 4ir/3pL 3 (15) 

where p and L are respectively mass density and the radius of the ECAL detector. We chose L ~ 
2OA511, where A511 is the attenuation length of 511 keV 7's. The relevant parameters for different types 
of materials used in ECAL's are listed below in Table 1. 

The required mass is minimal for a BGO ECAL due to its high effective Z (remember that the photo- 
absorption cross-section a ~ Z 5 ). Another important feature of BGO's is that they are not hygroscopic, 
thus, no additional dead material has to be introduced. For the crystal wrapped in aluminized mylar 
the light yield was measured to be 200±14 photoelectrons/1 MeV. This results in a probability of zero 
energy detection due to Poisson fluctuation of the number of photoelectrons, to be less than 10~ n for 
the zero energy signal defined as events with energy deposition less than 50 keV|2J]. A disadvantage in 
using BGO crystals is the 300 ns decay time of their scintillations, which may cause a drop of efficiency 
for the signal events due to the pile up effect. However for the source intensity planned to be used in the 
experiment this drop of efficiency is calculated to be about 15%. Thus, these results justify the selection 
of BGO as the 7-detector. The crystals, which have been lent to us by the Paul Scherrer Institute 

2 We would like to thank Dr.S.Asai for his ad vice 151. 
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ECAL type 


BGO 


Mat 


Csl( 11) 


Sc plastic/liquid 


Attenuation length, at 511 keV 


~1 cm 


~2.5 cm 


~1.9 cm 


~10 cm 


density, g/cm 3 


7.1 


3.6 


4.5 


1.0 


Required ECAL mass, kg 








~OOOlU 


Light yield, 7V 7 /511keV 


~ 4 • 10 3 


~ 20 • 10 3 


~ 10 • 10 3 


~ 2 • 10 3 


Hygroscopic 


no 


yes 


yes/slightly 


no 



Table 1 : Comparison between different types of ECAL 



(Villigcn, Switzerland), have a hexagonal shape with a length of 20 cm and an outer diameter of 5.5 cm, 
their original wrapping is a 0.75 mm thick teflon. In order to reduce this amount of dead material, the 
inner ring of BGO's has been wrapped in a 2 /im thick foil, aluminized from both side with 1000 A thick 
layers. The required number of crystals (~ 100), determined with the simulations, provides an almost 
isotropically uniform thickness of 20-22 cm of BGO. 

2.5 Monitoring 

Every crystal is equipped with a LED, fired with a pulser, to monitor the crystals response. The gas and 
light tight box is thermally stabilised to ±1°C. The HV for every PM is measured periodically during 
the run. 

2.6 Data Acquisition 

The data acquisition is started by the gate generated from the coincidence between the two XP 2020 
of the scintillating fiber. The gate has a length of 3.3 /is in order to suppress the probability of o-Ps 
decaying after this time to the level of 10 -11 . The energy of each crystal and the two XP 2020 signals 
from the fibers are recorded with four 32 channels CAEN QDC's v792. The timing between the signal 
from the fiber and the crystals of the endcaps and the inner most ring of the calorimeter is measured by a 
LeCroy 1182 TDC. The VME crate is interfaced to a PC with a National instruments VXI-MXI2 system 
and the data are acquired with Lab View. The monitoring is done by another computer using ROOT, 
not to slow down the data acquisition. 

3 Background 

The main background sources, summarized in Table 2, can be divided into three different types: 

• A real trigger accompanied by annihilation energy losses. 

— Hermiticity: as pointed out in the previous section (Table. 1), the 511 keV photons have an 
attenuation lenght of about 1cm in BGO. For 20 cm BGO the escape probability is about 

io- 9 . 
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— Absorption of photons in dead material: some energy of the photons can be absorbed either 
in the target or/and in the wrapping of the crystals or/and in the fiber. 

— Absorption of photons in the trigger energy window: the resolution of the endcap at 1.27 MeV 
is about 16% FWHM. If one (or more if e.g. one photon is backscattered in the target) of the 
positronium decay or direct annihilations photons overlap with the 1.27 MeV in a crystal it 
can be absorbed in the trigger energy window due to the resolution of the crystals. 

• A fake trigger, with no positron in the detector, can be produced by the EC 1.27 MeV photon. 
This gives two sources of background: 

— the 1.27 MeV photon can scatter with an electron in the fiber and the electron can deposit some 
fraction of its energy in the fiber and opens the gate for the DAQ. In this case the background 
level is about 1CP 6 . To reduce it, a trigger requires the coincidence between the signal from 
the fiber and the signal from the opposite (to the direction of the positron entering the fiber 
and the aerogel) endcap counter. This reduces the background by one order of magnitude. 
The next step is to suppress the background related to the 1.27 MeV photon scattering in the 
material, which surrounds the source, by depositing the 22 Na directly on the fiber. This lowers 
the background by another order of magnitude. The remaining background is associated with 
the 1.27 MeV EC photon scattering back to the trigger BGO (this photon will have between 
200-300 keV energy) and a Compton electron multiple scattering in the fiber hitting the same 
trigger counter. This seems to be the most dangerous background for our experiment. A 
possible solution, which is under study, will be discussed in Section 6. 

— The other possibility is that this photon may accidentally coincide with a fake positron trigger 
generated by the PM noise. 

• Physical backgrounds: 

— single photon or photonless decays of o-Ps ions: o — Ps~ — > e~ + 7, o — Ps~~ — > e~ + e~ 

— photonuclear absorption of annihilation photons accompanied either by photo-neutrons or by 
excitation of long lived nuclear states. 

The detailed evaluation of this background is in progress; some preliminary estimates show that it 
is less than 10~ 10 . The o-Ps produced in an excited state, immediately de-excites to the ground 
state[22] due to the high collisional (10 4 collision/lifetime) quenching in the S1O2 ■ Since the 
experiment is not performed in vacuum the positron wave function is assumed to overlap always 
with an electron wavefunction. The largest lifetime of the positrons in our experiment is assumed 
to be those for the three photon decay of o-Ps, namely 142 ns, thus making the probability of the 
positron not to annihilate during the 3.3 /zs gate negligible. 

4 Monte-Carlo 

The feasibility of the experiment has been studied with a Monte Carlo simulation based on the code 
Geant 3. The simulation was cross-checked and tuned with our previous search for an exotic three-body 
decay of o-Ps The ineflicieny for 511 keV 7 detection in the measurements was predicted to be 

(1.0 ± 0.2) x 10~ 4 . The comparison with the data shows an agreement within 30%. The low energy part 
of the spectra, shown in Fig. 4, is quite well reproduced, nevertheless more work has still to be done. The 
losses due to light transportation, self absorption and also the non uniformity of the crystals have to be 
included. 



5 Sensitivity 

As mentioned before, the trigger is given by the coincidence between the signals from the two PM of 
the fiber. Off-line, the 1.27 MeV photon is selected in the trigger counter inside the energy window 
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Figure 4: Comparison between MC simulation (points) and data of the experiment (dashed line) 



S 



Background source 


Expected level 


Hermiticity 
Dead Material 
Absorption in trigger energy window 


< io- iU 

2 x 1CT 10 
10" 9 


Compton EC photon 
Accidental noise and EC photon 


1.2 x 1CT 8 
3.2 x 1CT 11 


Physical backgrounds 


1CT U) 



Table 2: Expected background level 



1170-1470 keV (see Fig. 5). In order to reduce the accidental background (see previous section) the time 
between the coincidence and the 1.27 MeV 7 is required to be inside the time window as shown in Fig. 
5. A cut on the energy of the PM is applied to reduce the accidental noise. The distance between the 
source and the trigger BGO has been choosen to be 2 cm as a compromise between the trigger rate and 
the background from the Compton EC photon. The expected trigger rate for a source intensity of 50 
kBq is 

Trigger rate — — ■ Positrons emission rate ■ 5% • Fiber ef ficiency (16) 

and knowing the fraction of o-Ps in the spectra, which is about 5% 

Number of o — Ps/day = Trigger rate ■ o — Ps fraction ■ 86400s ~ 3 x 10 6 (17) 

Using Poisson statistics and assuming no background event observed, one can calculate the sensitivity of 
the experiment |23j after 60 days of running: 

Br(o - Ps -> invisible) = 2.3/1.8 x 10 8 ~ 1.3 x 10~ 8 (18) 



6 Conclusions 

With the presented type of calorimeter a background free limit for the Brio — Ps — » invisible) of ~ 10~ 8 
is reacheable in about two months of data taking. Although not all the possible background sources have 
been identified yet, the dominant background seems to be associated with the EC photon that scatters in 
the fiber back to the trigger BGO counter. The accompanying Compton electron can deposit ~ 50-100 
keV in the fiber, which is sufficient to trigger the coincidence of the PM's. Since the energy deposition in 
the trigger counter is almost equal to the total initial energy of the 1.27 MeV photon, it is quite difficult to 
reject this type of events with the given resolution of the crystal. One possible way to discriminate these 
e~7 pairs is to use a thin plastic scintillator in front of the trigger BGO, with a faster timing response 
than the one of the BGO's. This work is in progress. 
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